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The transient dynamics of long overlap Josephson junctions in the frame of the 
sine-Gordon model with a white noise source is investigated. The effect of noise 
delayed decay is observed for the case of overdamped sine-Gordon equation. It is 
shown that this noise induced effect, in the range of small noise intensities, vanishes 
for junctions lengths greater than several Josephson penetration length. 
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I. INTRODUCTION 



Josephson junctions are physical systems 
with nonlinear dynamics, which are inter- 
esting to investigate both from experimen- 
tal and theoretical points of view. This is 
also in view of numerous applications of su- 
perconductive devices based on Josephson 
junctions, such as RSFQ devices, qubits, 
SQUIDs, etc.. [Barone & Paterno, 1982; 
Likharev, 1986]. There are, in fact, a lot 
of open problems in Josephson junction's dy- 
namics, due to its nonlinear character. For 
some devices, as the RSFQ, minimization 
of the switching time is required for better 
performance [Pankratov & Spangolo, 2004]. 
From this point of view, it is interesting to in- 
vestigate the influence of thermal noise on the 
statistical properties of distributed Joseph- 
son junctions. It is also necessary to note 
that, currently, all Josephson junctions are 
manufactured with the use of optical and 
electron-beam lithography [Dorojevets, 2002; 
Makhlin et al, 2001], and can always be 
considered as distributed. Moreover macro- 
scopic quantum phenomena in long Joseph- 
son junctions have attracted a lot of exper- 



imental and theoretical work recently [Wei- 
des et al, 2006; Mertens et al, 2006; Alfi- 
mov & Popkov, 2006; Kim et al, 2006; Fistul 
et al, 2003]. These junctions are character- 
ized by one or more dimensions longer than 
the Josephson penetration length or depth 
[Barone & Paterno, 1982]. 

It was shown in different physical systems 
that thermal fluctuations can considerably 
increase the decay time of unstable states 
[Agudov & Malakhov, 1995; Malakhov & 
Pankratov, 1996; Agudov & Malakhov, 1999] 
and the lifetime of metastable states [Man- 
tegna & Spagnolo, 1996; Agudov & Spag- 
nolo, 2001; Spagnolo et al, 2004], producing 
a nonmonotonic behavior of these quantities 
as a function of the noise intensity. These are 
the effects of noise delayed decay (NDD) of 
unstable states and noise enhanced stability 
(NES) of metastable states. Both noise in- 
duced effects are due to the nonlinearity of 
the potential profile and are enhanced by the 
inverse probability current. The NDD effect, 
in particular, consists in the increase of the 
mean switching time (MST) of a Josephson 
junction due to the influence of noise. The 
main interest to analyze here long Josephson 



r Distributed bias current 




V I /'V 






XXVAWVA 







i — Superconductors 
FIG. 1 The structure of distributed Josephson 
junction of the "overlap" geometry. 



current density i is normalized to the crit- 
ical current density of the junction, if(x,t) 
is the fluctuational current density, T is the 
®" ^^ntifluxon normalized magnetic field, L = l/Xj is the di- 
mensionless length of the junction. Further, 
we will only consider the cases when i > 1 
and T = 0, so in the potential profile there 
are no metastable states at all. Also let us 
consider only the case of homogeneous bias 
current distribution, when % is constant along 
the junction. In the case where the fluctua- 
tions are treated as white Gaussian noise with 
zero mean, and the critical current density is 
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fixed, its correlation function is 



junctions (LJJ) is the existence of noise de- (i f (x,t)i f (x' ,t')) = 2j6{x-x')S(t-t'), (3) 
layed decay (NDD) effect for the restricted 

range of parameters, which is important for where 
practical applications. 



7 = It/(J c Xj) 



(4) 



II. THE MODEL 

In the frame of the resistive McCumber- 
Stewart model [Barone & Paterno, 1982; 
Likharev, 1986] the phase difference of the 
order parameter ip(x, t) of a long Josephson 
junction of the overlap geometry (see Fig. 1) 
is described by the sine-Gordon equation 



dt 2 dt dx 2 



i - sm((f) +ij(x,t)(l) 



with the following boundary conditions 
d<p(0,t) d<p(L,t) „ 



dx 



dx 



(2) 



Here the time and the space are normalized 
to the inverse characteristic frequency uj~ y 
and the Josephson penetration length Aj, re- 
spectively. The penetration length gives a 
measure of the distance in which dc Joseph- 
son currents are confined at the edges of 
the junction. (3 = 1/a 2 is the McCumber- 
Stewart parameter, a = uj p /uj c is the damp- 
ing, with C the junction capacitance and 
Up = (2eI c /hCy/ 2 the plasma frequency cor- 
responding to the critical current I c . The bias 



is the dimensionless noise intensity, J c is the 
critical current density of the junction, I T = 
2ekT/h is the thermal current, e is the elec- 
tron charge, H = h/2ir with h the Planck con- 
stant, k is the Boltzmann constant and T is 
the temperature. It is possible to obtain the 
expressions ([3]) and (j3j) through the renormal- 
ization of the general formula for the noise 
intensity, derived for the fixed total critical 
current I c [Castellano et al, 1996; Pankra- 
tov, 2002]. 

Initially, the whole phase "string" f(x, 0) 
is located at the position <p , so there are no 
vortices in the junction. The dynamics of 
such a phase "string" is similar to that of a 
real long string falling down on a tilted wash- 
board potential addressing the general prob- 
lem of the diffusion of an elastic string on a 
tilted periodic substrate [Cattuto & March- 
esoni, 1997]. For our case i > 1, such initial 
phase is considered to be located in the inflec- 
tion point of the potential profile (fo = tt/2. 
It is possible to prepare such a state by a fast 
switching from the starting superconductive 
state, corresponding to a metastable state, to 
the resistive one, which becomes then the ini- 
tial configuration of the system. Due to the 
unstable character of the initial position (see 
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FIG. 2 Qualitative illustration of the potential 
profile u(ip) = 1 — cos(yj) — itp, with i = 1.2, 
and phase string (p(x, t) in the process of escape 
from the initial position. 



Fig. 2), the phase will begin the determinis- 
tic motion along the potential profile. The 
MST is defined as the mean time of phase 
<p existence in the considered interval [— n, n] 
[Malakhov & Pankratov, 2002], where P(t) 
is the probability that the phase is located in 
the initial interval 

T T dp(t) 

t= J tw(t)dt = J P(t)dt, w(t) = -^-A 



Here w(t) is the probability density of the 
escape times. The probability Pit) is evalu- 
ated numerically in the following way: if at 
the given time t > the realization of <p(x, t) 
is within the interval [— tt, it], the probability 
for the corresponding realization is one, oth- 
erwise it is zero. After N realizations, the 
average, both over the realizations and the 
spatial coordinate x from to L, is taken 
and, finally, the required probability P(t) is 
obtained. 

Numerical solution of the sine-Gordon 
equation ([I]), with boundary conditions ([2]), 
is performed by using the implicit finite- 
difference scheme [Zhang, 1991], but taking 
care to insert correctly the noise intensity 
[Federov & Pankratov, 2007]. Typical val- 
ues of the simulation parameters are Ax = 
At = 0.1 — 0.02, for the spatial and time dis- 
cretization steps, and R = 10 3 — 10 4 for the 



number of realizations. 

III. NOISE DELAYED SWITCHING 

First, it is interesting to check the limit- 
ing transition of a long junction to a point 
one for small junction lengths L C 1, For 
a point junction, in the case of large damp- 
ing /? C 1, the MST was found analytically 
[Malakhov & Pankratov, 1996] for an arbi- 
trary value of the noise intensity. That for- 
mula, however, was obtained for the case of 
constant critical current I c , and the noise in- 
tensity in that case was: 7 S = It /la where 
I c = J C X jL. This means that we have to 
scale the noise intensity as 7 S = j/L, with 
7 given by Eq.flU). We get, therefore, the 
following closed expression for the MST in a 
long Josephson junction with constant criti- 
cal current density 

T=-\r e-fi'M f ef^dcpdx 

P(p2 POO ~\ 

+ / e siv)Lh d v - e- f ^ L ^d(p\, (5) 

J if l J ip2 J 

with 

f((p) = cos 99 + iip . (6) 

In Eq. (JSJ)) y?o is the coordinate of the ini- 
tial delta-shaped distribution, and (pi^ are 
the boundaries of the interval, delimiting the 
potential well or metastable state of the su- 
perconductive state, before the fast switch- 
ing. To prevent misunderstanding with the 
use of such a renormalization procedure, one 
should consider the general expression for 
the noise intensity for fixed critical current 
[Castellano et al, 1996; Pankratov, 2002] and 
substitute the required bias current density 
in it. In Figs. E| H the behavior of MST 
clS db function of the dimensionless noise in- 
tensity 7 and the junction length L, respec- 
tively, is shown for two values of the bias cur- 
rent. The agreement between the theoreti- 
cal results, obtained from Eq. (|3j), and the 
numerical simulations of Eq. (pQ), for a long 
junction, is very good not only in the limit- 
ing case L < 1, but even up to L ~ 1. In 
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Formula (5), 


= 1.2 


Formula (5), 


= 1.5 
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FIG. 3 MST versus the noise intensity 7. 
Comparison between the theoretical results 
of Eq. ([5]) and the numerical simulations of 
Eq. (JTJ), for two values of the bias current, 
namely i = 1.2, 1.5. Here the dimensionless 
junction length is L = 0.1. 

Figs. El M the semilog plot of the MST versus 
the noise intensity 7, for different junction 
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FIG. 4 MST versus the dimensionless junction 
length L = l/Xj. Comparison between the the- 
oretical results of Eq. ([5]) and the numerical 
simulations of Eq. (JT]), for the same values of 
the bias current i of Fig. 3. Here the dimen- 
sionless noise intensity is 7 = 0.005. 

lengths, is shown. The NDD effect is present 
for junction lengths up to L > 5, while for 
greater lengths it completely vanishes. This 
peculiarity is very important for RSFQ de- 
vices. This disappearance of NDD is due to 
the effective decrease of the noise intensity in 
comparison with the point junction. 
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FIG. 5 Semilog plot of MST versus 7, for 
a long Josephson junction with i = 1.2 
and different junctions lengths, namely L = 
0.1,0.5,2,10. 
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FIG. 6 Semilog plot of MST versus 7, for a 
long Josephson junction with i = 1.5 and the 
same different junctions lengths of Fig. 5. 

For the same temperature level, the effec- 
tive noise intensity 7 S = 7/L for the long 
Josephson junction will be smaller by a fac- 
tor 1/L in comparison with the point junc- 
tion. So, for LJJ with rather large lengths 
(L > 5) the noise intensity will get out of the 
area of the NDD effect. 

For example, from our simulations, we 
find that for great values of the noise inten- 
sity 7 ~ 100, the NDD will exists even for 
LJJ with lengths L ~ 10, but such a range 
of 7 is not interesting from practical point 
of view. In other words, for large lengths 
of the distributed Josephson junctions, the 
random force Ft ~ -Jj becomes negligi- 



4 



bly small in comparison with the determin- 
istic force, caused by the potential profile 
u(ip) = 1 — cos(v?) — itp. 

The standard deviation (SD) a of the 
switching time is defined as 

^(i 2 ) -r 2 , (t 2 ) = rt 2 w(t)dt. 



a 

For small noise intensities we can obtain an 
expression of a((p ) for a long Josephson junc- 
tion, by using the noise intensity renormal- 
ization, as before in Eq. and the asymp- 
totic expression derived for a point Josephson 
junction [Pankratov & Spagnolo, 2004] 

*M = -y/&Y/L) [F(^ ) + /3^o)] + ..,(7) 

UJ r . 



h{x) 
h{<Pa) 



2 . ( ita.n(x/2)-V 

— - — ■■ arctan ( 



fV2 



l/(2(i- 

cos(x)/i(x) 
(sin(x)— i) 3 



sin x) 2 ), 



dx. 



2(sm(z)-i) 3 _ 

In the following Fig. 7 the SD of the switching 
time versus the junction length L is shown, 
for two values of the bias current, namely 
i = 1.2, 1.5. The theoretical results of Eq. (f7|) 
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FIG. 8 Standard deviation versus the dimen- 
sionless noise intensity 7 for four junction 
length, namely L = 0.1,0.5,2,10. Here the 
bias current is i = 1.2. 
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FIG. 7 SD versus the dimensionless junction 
length L for two values of the bias current, 
namely i = 1.2,1.5. The theoretical results of 
Eq. ([7]) are compared with the numerical sim- 
ulations of Eq. dl]). Here the noise intensity is 
7 = 0.01. 

show a quite good agreement with the nu- 
merical simulations of Eq. ([1]). In Figs. 8, 9 



FIG. 9 Standard deviation versus the dimen- 
sionless noise intensity 7 for the same junc- 
tion lengths of Fig. 8. Here the bias current is 
% = 1.5. 

the numerical simulations of the SD versus 
the noise intensity from Eq. 01) , for two val- 
ues of the bias current and for different junc- 
tion lengths, are shown. All the behaviors 
of Figs. 7, 8, 9 are evaluated in a parameter 
range which is characteristic for RSFQ de- 
vices. The SD decreases with the length of 
the junction. This is because the dynamics 
controlling the switching event goes from a 
noise-induced regime, at very short junction 
lengths, to a deterministic regime, caused by 
the decrease of the effective noise intensity 
for very long junctions. The noisy regime is 
well visible in Figs. 8, 9 for the case L = 0.1, 
for which the a increases with the noise in- 
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tensity. The deterministic regime, dominated 
by the potential profile, produces the asymp- 
totic constant values of the SD (the curves for 
L = 0.5, 2, 10 in Figs. 8, 9). We note that for 
a higher bias current (Fig. ) the SD reaches 
lower asymptotic values, as we expect. 



IV. CONCLUSIONS 

In this work a study of the transient dy- 
namics of long Josephson junctions, with bias 
currents greater than the critical current, is 
presented. We have found that, for the case 
of constant critical current density, the mean 
switching time has a maximum for small 
noise intensities and junction lengths. This 
is the noise delayed decay effect. This effect 
disappears for junctions with dimensionless 
lengths L > 5, which can be important for 
the design of RSFQ devices based on Joseph- 
son junctions. The disappearance of NDD 
effect is explained by the decrease of the ef- 
fective noise intensity with the length L of 
the junction, by considering the scaling fac- 
tor 1/L. As a consequence the system gets 
out from the range of 7 values suitable to 
observe the noise delayed decay effect. Fi- 
nally we observe for the MST and its SD, in 
the range of small noise intensities, a good 
agreement between the theoretical behaviors 
(Eqs. ([S]) and ([7])) and those obtained from 
numerical simulations of Eq. ([T]). 
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